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We study grand-canonical and canonical properties of the model of 
branched polymers proposed in We show that the model has a fourth 
order phase transition and calculate critical exponents. At the transition 
the exponent 7 of the grand-canonical ensemble, analogous to the string 
susceptibility exponent of surface models, 7 ~ 0.3237525... is the first 
known example of positive 7 which is not of the form 1/n, n = 2, 3, . . .. 
We show that a slight modification of the model produces a continuos 
spectrum of 7's in the range (0, 1/2] and changes the order of the transi- 
tion. 



Introduction 



The problem of summing over random geometry appears in many areas of 
modern physics such as string theory, quantum gravity, membranes and oth- 
ers. The problem is known for being hard to study both analytically and 
numerically. Therefore one frequently uses simplified models. In many cases 
branched polymers capture some essential futures of more complicated mod- 
els. Being simple and solvable they offer us insight into such issues like value 
of the string susceptibility exponent, correlation functions, renormalization 
group, sum over genera etc 

In the present paper we study a model with a coupling to the branching of the 
polymers The model has two regimes. For small values of the coupling the 
exponent 7=1/2 while for larger values it is negative and depends linearly 
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on the coupling. We extend the analysis from the original paper localize 
the transition and find that at this point 7 ~ 0.3237525.... 

We also study the canonical ensemble and show that the model exhibits a 
fourth order phase transition. A slight modification of the model allowed us 
to obtain a continuos spectrum of 7's between and 1/2. In the canonical 
ensemble this generalized model can exhibit a phase transition of an arbitrary 
order equal or higher than 2. 



The model 

We consider an ensemble of planar rooted trees generated by the recursion 
relation shown in Fig.Q. Each bubble in the figure denotes the grand canon- 




Fig. 1. Partition function 

ical partition function Z. By iterating the relation one produces all possible 
planar rooted trees. Each tree contributes to the partition function the factor 
e~'^", where n is the number of vertices in the tree. Each vertex additionally 
contributes a weight factor depending on its degree, ie the number of branches 
which meet at vertex. In Fig.|I| we denoted these weight factors by t^, where k 
is degree of vertex. In this paper we analyze the model with tk = k~^, where 
/3 is an external parameter. The equation from Fig. |l] can be written as : 

which generates following partition function (in the following we drop the 
arguments of the Z function) 

Z=Y^ e-/.n(T)g-/3i?(T) (2) 

where T denotes the ensemble of all trees. The function E{T) plays the role 
of energy which the tree contributes to the ensemble and is equal : 
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E{T) = J2^nk{v) (3) 

The sum goes over all vertices of the tree T, and k{v) is degree of vertex v. 
By summing first over trees with fixed number of vertices we rewrite (^) as 



E 

n=l 



-fin 



n=l 



(4) 



The coefficients z„ of the series at e"^"' depend on (3 and correspond to a sum 
of e~^^^'^^ over the ensemble of all trees with n vertices. They are canonical 
partition functions for the ensembles of trees with n nodes. We have introduced 
convenient variables which we will use later on : free energy density fniP) for 
the system with n vertices defined from the expansion coefficients : Zn = e"'-^" 
and thermodynamic potential which we call free enthalpy (p = —InZ. 

For large positive /3 the canonical partition function is dominated by branched 
polymers with minimal energy. Such branched polymers maximize the number 
of vertices in the last generation ie vertices with degree k = 1. They look like 
short bushes. In the extreme case such bushes have one branched vertex with 
n branches while all other vertices belong to the last generation. The energy of 
such branched polymers is In ra. When (3 becomes smaller the entropy becomes 
important and more and more elongated structures appear in the ensemble. 
At certain /? the system undergoes a fourth order phase transition associated 
with the change of the regimes. Finally if one allowed for negative f3 then 
for large negative f3, trees with large energy would mainly contribute to the 
canonical ensemble. The maximal energy E for the system with n nodes is 
n In 2 and comes from the chain structure. 



Exponent 7 



The series representing the grand canonical partition function Z is regular 
at e~'^ = 0. The behaviour of Zn in the thermodynamic limit n ^ 00 can be 
deduced from the behaviour of the function Z at the radius of convergence 
e~'^°, when fi ^ fiQ, which corresponds to the first singularity of Z encountered 
when fi comes from 00. At this singularity free enthalpy as a function of 
chemical potential behaves for a given (3 as : 



- 00 ~ (/i - l^oY ^ (5) 

where 0o, l^o and 7 depend on [3. The parameter /iq is the critical value of 
the chemical potential and 7 is the universal exponent defining a type of 
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singularity which accounts for the universahty class a model belongs to. Its 
counterpart in the surface models is called the string susceptibility exponent. 
As will be seen later 7 is equal or less than 1/2 and therefore the function on 
the right hand side goes to when /i approaches /iq. The value of free enthalpy 
aX fi = ^0 is denoted by 0o = 0(/io)- The singularity (|^) corresponds to the 
following asymptotic behaviour of the coefficients Zn in the limit of large n : 



ZniP) ~ n^-'e^>'" (6) 
In the thermodynamic limit n ^ 00 the free energy density is foo{f3) = fJ'oiP)- 
The series Z{fi) is given by an indirect relation (||) which can be rewritten 

/i(0,/3) = lnL^(e-^) + 20 (7) 

where Lplz) = Z^J^Li -2"/'^'^ is a special function whose basic properties, fre- 
quently used in this paper, are listed in the appendix. Throughout the paper 
we use convention for functions which have (3 as second argument that prime 
denotes derivative at fixed j3. The function 4>{n) is singular at fiQ when either 
/^'(0o) = or /i(0o) itself is singular at (pQ = (t>{fio)- The first condition reads 



/-m/3) = ^^ = (8) 



which, by virtue of (p2D, simplifies to : 



2L^(e-'^«) =L^_i(e-<^«) (9) 

The solution of @ is shown in Fig. ^ The function 0o approaches its critical 
value (pc = when P [3c- At this point the argument of Lp : e~^° becomes 
unity and is equal to the radius of convergence of L^. For j3 > Pc the critical 
value of (j) stays at 0c = and corresponds directly to the singularity of the 
map /3) . Inserting e~'^° = 1 to the (H) we get the critical value of /3 : 



2C(/3c) = C(/5c - 1) (10) 
which can be found numerically Pc = 2.4787508.... 

Inserting this back into /i we get the critical value of the chemical potential at 
the radius of convergence of the series (^) which also corresponds to the free 
energy density in the thermodynamic limit foo{P) '■ 



4 




Fig. 2. Critical free enthalpy 



/oo(/3) = /io = /i(0o(/3),/3) 



lnL^(e-<^«(^)) + 2<^o(/?), 13 < I3c 



lnL^(l 



lnC(/3), 



/3> A 



We can now compute the exponent 7. We separately consider three cases 
(3 < (3c, (3 = (3c and (3 > (3c- We invert the function with respect to 

0, for near to 0o- For P < Pc have : 



and from there 



+ . . . 



(12) 



9 \ 1/2 

4) (/i-/^o)^/^+... 



The second derivative /Xq is finite in this range. For (3 
starts directly from the singular term : 



(13) 

Pc the expansion 



fi-fio = w{Pc){(pf^-' + ... 



(14) 



because fi' = at this point and Pc is between 2 and 3. The coefficient 
w{P) = r(l — P)/({P). Inverting this with respect to (p we obtain : 



cf> = w{PcY/^^--'\f^-f^oY^^^--'^ + 



(15) 



Finally, for P > Pc the expansion of fi around /xq starts with a nonvanishing 
linear part : 
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II- fio = Aio0 + . . . + w{f3)(f)'^ ^ + . . . 



(16) 



The dots between the terms (p and 0^ ^ correspond to integer powers 0" with 
n less than /5 (see (^4|)). Inverted with respect to the equation yields : 



— (/i-/io) + ... 



Comparing the singularities with 



ifi - fiof-^ + ... 
we obtain : 



(17) 
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1 

2' 

f3c-' 
2-/?, 



0.3237525..., = 



c 



c 



Canonical ensemble 



Properties of the canonical ensemble in the thermodynamic limit are described 
by the free energy density foo{P)- It is a continuos function of /5 but it has 
singularity at /5 = /3c. To extract the dominating singularity of free energy 
for Ap = {Pc — A) — * 0+ let us first consider the critical free enthalpy 0o(A) 
defined by the equation (||). Using (0) we can for small (3 and 0o expand the 
equation (P) : 



~aAf3 + b<f)^o^-^ + ... 

up to o(A/3) and o(0(^^"^ 
('(pc - 1) = 3.4950..., b 
equation gives 







(19) 



The expansion coefficients are a = 2('{Pc) ^ 
-r(2 - pc) = 3.5548.... Solved for 0o the last 



iP) = + 



(20) 



in the limit Ap 0+. The exponent r] = 3 - l/{Pc - 2) = 0.911231..., and 
the coefficient c = (a/fe)^/(3-r,) ^ 0.9666.... We can now insert 0o into ([IT]) and 
compute singularity of foc{P)- The first derivative is : 



L^(e-'^o) 



(21) 
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The term containing dfi/d(f) vanishes identically for 0o ®- Now we can ex- 
pand the last expression in small 00; and then in small A(3. The dominating 
singularity for Af3 O"*" comes from the first nonvanishing term containing 
00- This term is linear in 0o and introduces the singularity A/3^~^ : 



fLW) = ^-^ + dA(3'-^ + ... (22) 



The coefficient in front of the singular term : 




((lie) \ ai3c) 

is obtained by expanding (^T]) and using (|34D . The first term on the right hand 
side of (^) matches the function for [3 > [3c (JH]) and solutions come smoothly 
one into the other. The contribution of the singular term (A/5)^~'' vanishes at 
A/3 = as well as for the second and the third derivative of the free energy 
density f^o- For the fourth derivative the singular term blows up at [3c '■ 



'^m-'^m-e^r^,... (23) 

where e = (3 — ?7)(2 — r7)(l — ifjd = 0.5045.... This shows that the model has 
a fourth order phase at Pc in the canonical sector with a singularity given by 
the exponent r]. 



Non— universality 



The dependence of the exponent 7 on the critical value of the coupling suggests 
a non-universality. By changing slightly the model one should be able to 
change the value of Pc without spoiling the formulas ([T8|) . In this way one 
can change the value of the exponent jc- We will show below that indeed a 
simple modification of the model produces a whole continuos spectrum of 7's. 

We modify the model by letting the weight ti, coupled to the number of points 
of the last generations, vary in contrast to the model defined by (|I|) where it 
was set to 1. Now ti can take arbitrary value greater than zero. The value of 
ti affects the number of points in the last generation. Decreasing/increasing 
ti for a given P favours/suppresses polymers with many points in the last 
generation. By changing ti we can shift the critical value of P where the phase 
transition takes place. On the other hand this modification does not change 
the singularity of the function on right hand-side of ([^) . 
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After introducing the modification of the weight ti the equation (0) becomes 



(3) = ln{L^(e-'^) + {t, - 1) e"-^} + 



(24) 



Repeating the analysis from the previous sections we first find the critical line 
00 as a solution of /i' = : 



(l-ti)e- 



(25) 



for P less than Pc and for /? greater than Pc- The position of the critical 
point is given by : 



2C(/5c) - CWc - 1) = 1 - ti 



(26) 



The value Pc as a function of ti is shown in Fig.^. For ti +oo, the value 
of Pc goes to 2. For ti ^ 0, Pc goes to infinity. As first, consider the case 




Fig. 3. Phase diagram of the extended model 

2 < Pc < 3 corresponding to t > = 1 + C(2) - 2C(3) = 0.24081 . . . . 
Repeating the same analysis as for the original model we get that 7 is given 
exactly by the same formula ([l^) with Pc replaced now by the current value 
Pciti)- At the critical point we have •jc = {Pc ~ '^)/{l3c — 1) which changes 
contiuosly between 0, 1/2 when Pc goes from 2 to 3. When Pc exceeds 3, 
which is for ti < t*, the situation changes. The expansion ([T^) starts then not 
from the singular term with the power Pc — 1 but from the term /igV^ which 
after inverting gives the singularity ~ (/x — /io)^^^ and hence 7 = 1/2. 

To find the type of the singularity in the canonical ensemble we determine the 
behaviour of the critical free enthalpy for A/3 = {P — Pc) — 0^ . As before 
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first consider the case 2 < /^c < 3. Expanding the equation /i' = for small 
A/3 and small we get exactly the same equation as ([T9|) with exactly the 
same coefficients but calculated now at different (3c- The function 0o has then 
the singularity 0o ~ (A/3)^/'^^c'^^) whose power depends through Pc on the 
parameter ti and varies continuosly. 

The situation changes for /3c > 3 (ti < t^) because then additionally regular 
terms in 0o enter the expansion. They introduce powers of 0o smaller than 



~aAf3 + g<Po + ... + b4^-^ + ... = (27) 

where g = —2Q{l3c — 1) + C(/^c — 2), and a, 6 are as in the equation (|T^. 
Inverting the last equation for 0o as a function of A/3 we get : 



00 = -A/? + ... - ^"-^0-2) (28) 
9 9 9 

Inserting the singularity of 0o into analogously as we did in (|^) we obtain 
the singular part of the first derivative of the free energy : 



, ' (29) 

One sees that the order of the transition grows arbitraryly and the transitions 
gets softer when Pc approaches 2, corresponding to ti 0, or when Pc, 
tl —>■ oo. The strongest transition takes place for ti = tl for which /3c = 3 and 
the second derivative of the free energy diverges logarithmically ln(A/3). 



Discussion 



The model described above exhibits a manifold behaviour summarized on 
the phase diagram (ti,/3) (Fig.|]). The model has two phases with bush and 
tree like branched polymers. They are separated by a critical line on which the 
order of the transition and the critical exponents vary. It has a stable phase for 
P < Pc(ti) where the exponent 7 = 1/2 and it is not affected by a perturbation 
of couplings. This is a generic branched polymer phase and has realizations 
in surface models and higher simplicial gravity At the critical line Pc'{ti) 
the exponent 7 is unstable against small changes of the model parameters and 
reveals a continuos spectrum in the interval (0, 1/2]. 
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It is known that for c > 1 exponent 7 is not universal which means that 7 
depends not only on c but also on additional parameters controlling a specific 
model realization. For example, for gaussian field on random surface the ex- 
ponent 7 chang function on the number of fields c and the integration 
measure parameter ||^. A classification for a certain subclass of models was 
listed by Durhuus who showed that positive 7 is of the form 1/n n = 2,3, ... 
and in particular 7 = 1/3 is realized by multispin models. Our results do not 
prove that surface models with arbitrary 7 in the range (0, 1/2] do exist but 
surely make it more plausible. 

We believe that models of branched polymers are worth studying not only 
in the above mentioned context of universality of random surfaces but also 
as solvable models for fiuctuating geometry. Our present understanding of 
critical phenomena is mainly based on the spin type models on regular fixed 
geometry. One wonders if this intuition can be translated directly to fiuctu- 
ating geometry. In particular, there is a difficulty with defining correlation 
functions on fiuctuating geometry since the distance between points fiuctu- 
ates 0J^. A study of models as the one discussed here, especially around the 
phase transition, can hopefully give us a more profound understanding of those 
problems. 
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Appendix 

We shortly summarize properties of the function defined by the series : 

L,iz) = E (30) 

n=l 

For z = 1 and /? > 1 the function reduces to the Riemann Zeta function : 

= m (31) 
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Calculating derivative with respect to the argument z we get within the radius 
of convergence 1^1 < 1 : 



z"-^ = L,.,K-A (32) 

which also holds for z = 1 when j3 > 2. The function Lp{z) has singularity 
when 2; — > 1 of the type : 



sing L(,{z) ~ r(l -m- zf-\ /3 ^ 1, 2, . . . (33) 
where F is the Euler gamma function. 
Put together, the last equations give for x — : 



Lp{l -x) = m + ... + (-)"C(/3 - + r(l - ^)x^~' (34) 

up to o{x^^^). n is the largest integer less than (3 — 1. For integer /3's the 
dominating singularity is of the type x^~^ \n.x. 
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